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HEAT TRACE ASYMPTOTICS OF 
SUBORDINATE BROWNIAN MOTION 
IN EUCLIDEAN SPACE 

M.A. FAHRENWALDT 


Abstract. We derive the heat trace asymptotics of the generator of subor¬ 
dinate Brownian motion on Euclidean space for a class of Laplace exponents. 
The terms in the asymptotic expansion can be computed to arbitrary order and 
depend both on the geometry of Euclidean space and the short-time behaviour 
of the process. If the Blumenthal-Getoor index of the process is rational, then 
the asymptotics may contain logarithmic terms. The key assumption is the 
existence of a suitable density for the Levy measure of the subordinator. The 
analysis is highly explicit. 


1. Introduction 

This paper explores the correspondence between stochastic processes on R” and 
analytical objects: we pick a subordinator X t . i.e. an increasing process with inde¬ 
pendent and homogeneous increments, and let Bx t be the subordinate Brownian 
motion. It is obtained from a standard Brownian motion by introducing a local 
time given by the subordinator. Let A be the generator of the corresponding semi¬ 
group. We show that in the heat trace asymptotics of A we recover both geometric 
information about R" and probabilistic information about Bx, and X t . This is 
illustrated schematically in the following diagram and will be made precise in the 
next section. 


Stochastics 
Process on R” 
Bx t 


t~ x !P sup \Bx a \ 4' 

S<£ 


0, /3 > 2a 

oo, j3 <2a 


as t —>■ (P with Blumenthal-Getoor index 2 a 


Analysis 

Heat kernel 
„-At 


e ct J2 c k (a)t( n ~ k '>/ 2a 

-ATR(e~ At )~ k =° 

' ' OO 

-e ct c k (a)t k \ogt 


k =1 


heat trace asymptotics for t O - ^ - , 
log terms only for o: rational 


Three features differentiate our results from the literature. First, we consider 
processes that live in Euclidean space and not in compact domains. Second, we 
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observe that the heat trace asymptotics strongly depend on the (ir)rationality of a 
parameter that governs the short-term behaviour of the process. Third, our method 
allows to compute the asymptotics to arbitrary order. 

We consider a class of subordinators that is small enough to allow a fully tractable 
analysis yet large enough to show interesting behaviour. Roughly speaking, we 
work with subordinators whose Laplace exponent in its Levy-Khintchin form has 
a density with suitable asymptotic expansion near the origin and is of rapid decay 
at infinity. This class contains the relativistic stable process, which is important 
in applications in financial mathematics or quantum physics. We present this as a 
fully worked example in dimensions 2 and 3 demonstrating that the latter situation 
leads to logarithmic terms in the heat trace asymptotics. 

The essential technical features of our approach are the use of a global calculus 
of classical pseudodifferential operators on R n and a generalized trace functional 
on this algebra. The trace functional allows us to go beyond compact domains and 
consider processes with values in the whole of R". 

More precisely, for our class of subordinators the suitably shifted generator A = 
A — cl belongs to the algebra of classical S'G-pseudodifferential operators on R n 
m m and so do its complex powers A z and the heat operator e At . Using 
a generalized trace functional TR analogous to [ 23 ] on a suitable subalgebra we 
explicitly compute the regularized zeta function £ = TR(A~ Z ) and the generalized 
heat trace TR(e~ At ). The asymptotics of the latter can be expressed in terms of 
the pole structure of the regularized zeta function. The trace functional TR is the 
analogue of the Kontsevich-Vishik trace m m on closed manifolds which is used 
to study determinants of elliptic operators. These play an important role in many 
areas of mathematics and quantum physics, cf. [3Tj for comprehensive references. 

In a broader perspective, we argue similarly to noncommutative geometry 131: 
given a stochastic process we try to infer probabilistic information from the spec¬ 
trum of a strategically associated operator. 

Closely related to this investigation are [7 [26] that compute several terms in the 
heat trace expansion for the relativistic a-stable processes on a compact domain in 
Euclidean space. We comment on this in Section [5] and also show that our results 
agree with [7] and [26], cf. Remark |5. II The case of subordinate Brownian motion 
on closed manifolds is covered in mmmm- 

The use of pseudodifferential operators to understand Feller processes is the 
theme of the comprehensive series m [T9] 20] . Building on the seminal work of 
Hoh-Jacob-Schilling, the author describes a symbolic calculus for a class of pseu¬ 
dodifferential operators that appear naturally as generators of certain Markov pro¬ 
cesses. These operators have a more intricate structure than the pseudodifferential 
operators used in index theory yet allow for parametrices and the Fredholm prop¬ 
erty. Also, the Ruzhansky-Turunen theory of pseudodifferential operators on Lie 
groups has found successful application in Markov processes as developed by Ap- 
plebaum [3]. Either calculus appears, however, not suitable for our purposes due 
to the lack of a homogeneous symbol expansion. 

This paper is organized as follows. The following section recalls the basic nota¬ 
tion for symbols of pseudodifferential operators. Section [3] contains the statements 
of our key results with proofs given in Section [4] Section [5] contains a fully worked 
example also illustrating several calculations from the proofs. 
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2. Symbol spaces 


We briefly recall the symbol spaces of the theory of S'G-pseudodifferential op¬ 
erators on R", the reader is referred to HU ns Ea mu ng for a more detailed 
treatment. For the notation to be consistent in this paper we need to adjust some 
of the standard notation in pseudodifferential operators. 

We define the smooth weight function ( x) = (1 + |a;| 2 ) 1//2 and introduce the 
partial differential operators dj = . For a multi-index /3 = {f3\ £ INq 

we set = <9f 1 • • • d@ n . 

Denote by C' 00 (R") the set of smooth functions R" -A- (D and let <S(R") be the 
Schwartz space of rapidly decaying functions, i.e. u £ C°°( R n ) such that 

sup \x 1 d^u(x)\ < oo 

rc6]R" 

for all /3,7 £ Nq. Denote by 5'(R n ) the topological dual of the Schwartz space, 
called temperate distributions. 

As usual, let L 2 (R") be the Hilbert space L 2 (R") ={«£ ^'(R’^Kii, u) < oo} 
with inner product (u,v) = f u(x)v(x)dx and norm |M| 2 2 ( R n) = (u,u). 

Define the Fourier transform u(£) of a function u £ <S(R") by 


where rfx = (27r) n ! 2 dx and • denotes the standard inner product in R". The 
Fourier transform yields an isomorphism of 5(R") extending to an isomorphism of 
«S'(R n ). 

We now introduce the symbol class S r ’ r (R n ). 


Definition 2.1 ([233, Definition 2.1). The space S r,r (R ra ) of symbols of order 
(r, r') £ (D 2 is the set of smooth functions a : R n x R n —> (D such that for all 
multi-indices /3 ,7 £ Ng there is a constant Cp n with 


d^d2a(x,£) 


< G^, T (x) Rer '- | 7 l (e ) Rer " l/31 


for any (a;,0 £ R n x R n . We set S-°°’ r '( R n ) = n r . eR S' r '- r ''(R"). 


The pseudodifferential operator A with symbol a £ S r,r (R") is given by 

(1) Au(x ) = J e lx '^a{a:,£)u(£)dx 


mapping S{ R n ) -A- <S(R n ) and S( R n )' S( R")'. 

For symbols independent of x we recall symbol expansions and classicality of 
pseudodifferential operators. The classical symbols have asymptotic expansions in 
£ into homogeneous terms, we refer to Definition 2.2 of [23j for the general situation. 

Definition 2.2. Let <r(£) £ 5 ,T ' , 0 (R”) be a symbol independent of x. 

(i) We say that it has the asymptotic expansion 

OO 

~ ^Vr-kiZ) 

k =0 

where cr r -k £ G 00 (R" \ {0}) if the following holds. For a cutoff function 
X £ G°°(R") equal to 1 outside |£| > 1 and equal to zero for |£| < 1/2 we 
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have 

N—l 

^-E^KJeS^'r) 

fc =0 

for any N > 1. 

(ii) The symbol is classical if the a r -k are homogeneous in £, i.e. if oy_fc(r£) = 
T r_fc ay_fc(£) for r > 0. The set of classical symbols is denoted by 5^°(1R"). 

We also recall the notion of ellipticity and parameter-dependent ellipticity in the 
context of symbols that are independent of x. We refer to [53], Definition 2.3 and 
[24] , A.3, respectively, for the full generality of the definitions. 

Definition 2.3. Let <r(£) G S r,0 (lR n ) be a symbol independent of x with order 
r > 0. 

(i) We say that a is elliptic if there is a constant C > 0 such that 

C(£) r < *(£) 

for all (eE". 

(ii) Fix 0 < 9 < tt and let A = {re llp \r > 0 and 9 < ip < 2n — 9} be a sector in the 
complex plane. We say that a is A -elliptic if there is a constant C > 0 such 
that 

(a) er(£) does not take values in A for all £ G IR", and 

(b) |(A — cr(£)) _1 | < C{£,)~ r for A G A and £ G R". 

3. Statement of the key results 
W e state and motivate the assumptions and formulate the key results. 

Probabilistic prelude, cf. m mi- Let B t be a Brownian motion on ]R ra 
with characteristic function 3E(e l ^' St ) = e -t ^l for £ G R n and t > 0. Let X t be 
a subordinator on [0, oo) independent of i.e. an increasing Levy process with 
values in [0, oo) and Xo = 0 almost surely. The distribution of B\ t can be described 
in terms of the characteristic function lE(e^' Bx *) = 1 for t > 0, where / 

is the Laplace exponent in the probabilist’s convention or the Bernstein function 
in the analyst’s vocabulary. Also the generating function of the subordinator is 
E ( e -AY t ) = g—t/(A) for A > o 

Recall (cf. [5D], Definition 3.1) that a function / : (0,oo) —> R is a Bernstein 
function if / is smooth, /(A) > 0 and (—l) fc_1 /^ fc ^(A) > 0 for k G IN. Any Bernstein 
function can be represented in Levy-Khintchin form as 

nOO 

(2) /(A) = a + bX + / (l - e~ xt ) p(dt), 

Jo 

for constants a, b > 0 and p, a measure on (0, oo) such that / 0 °° t A 1 p(dt) < oo. 
The Levy characteristic triplet (a, b , p) uniquely determines /. 

We consider Bernstein functions whose Levy measure has a locally integrable 
density m with respect to Lebesgue measure, and we call this the Levy density. 
Moreover, we restrict ourselves to Bernstein functions of the form <[2|) with a = b = 0 
but indicate in Section [4] how the general case differs. 

One can canonically associate a semigroup T t with the process B\ t by defining 
[T t u](x) = lE x (u(Bx t )) acting on Schwartz functions u G S(R"). The generator of 
this semigroup is defined as the operator Au = lim^o TtU t ~ t with domain the set 
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of functions where this limit exists. The link with the probabilistic picture is that 
the generator acts as the integral operator 

(3) Au{x) = -(2ti r)-"/ 2 / (|£| 2 ) u{^ 

JR" 

where u £ <S(1R") and u denotes the Fourier transform of u. So if / is smooth and 
its derivatives decay sufficiently fast, then A is a pseudodifferential operator with 
symbol — /(|£| 2 ). We sometimes call this the symbol of the subordinate process. 
We recall the definition of asymptotic expansions of real-valued functions. 


Definition 3.1. Suppose that g : (0,oo) -A 1R is a function. We say that g(t) ~ 
SfcLo Pkt ak as t —>■ 0 + if pfe £ 1R, a/c t oo and 


lim t~ aN 
t-> o+ 


N 


g(t) - Y,Pk tak 


fc =0 


= 0 


for every N > 0. Analogously for t —> oo. 


The key assumption in this paper is the existence of a suitable Levy density. 

Hypothesis 3.2. Let /(A) = / 0 °° (l — e -At ) m(t)dt be a Bernstein function with 
locally integrable density to : (0, oo) —> 1R that has the following properties. 

(i) There is an a £ (0,1) such that m has the asymptotic expansion 

OO 

m(t) ~ t-^y^ 
k =o 

as t —> 0 + . 

(ii) to is of rapid decay at oo, i.e. TO.(t)t^ is bounded a.e. for t > 1 for all /3 £ 1R. 

(iii) to(0, oo) < 0 where to(0, oo) = / 0 °° ( m(t ) — po^ _1_ “) dt. 

Assumption (i) yields an asymptotic expansion of / for large A and assumption 
(ii) makes / smooth at the origin as it implies that f 0 t l m(t)dt < oo for any l £ N. 
Moreover, (iii) is a technical condition ensuring A — m(0, oo) with A as in ([3]) is an 
invertible and classical pseudodifferential operator. 

Example 3.3. From [29] we pick five examples of Bernstein functions that sat¬ 
isfy Hypothesis 13.21 The asymptotics of the Levy densities can be obtained using 
Taylor’s theorem. In each case, a £ (0,1) and c > 0. 

(i) /(A) = (A + 1)“ — 1. This has Levy density m{t) = r^i-a) and the 

asymptotic expansion 

™(t) - rH(! - * + + ■ ■') 

as t —> 0+. This Bernstein function describes the relativistic a-stable Levy 
processes, which is related to the relativistic Hamiltonian in physics, cf. [I] and 
to the Normal Inverse Gaussian distribution used in financial mathematics, 
cf. [8]. All key results are illustrated for the case a = 1/2 in Section [5] 

(ii) /(A) = A/(A + c)“ with m(f) = sln ( n7r T( 1 ~ a ) e -at^a- 2 ^ _|_ \ _ a ) anc [ asymp¬ 
totics 

m(t) ~ ^((i-o>)r(A) t -i_ a ' ^ + c(1 _ a , )t + c 2 (i a /2 _ !) t 2 + ...) 
where a' = 1 — a. 
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- 2^»/=» - ' Then 


(iii) /(A) = A ^1 — e 2 \A+ C ^ />/A + c and 777 (f) = - 

m(t) ~ 2 ^ rl ~ 1/2 (l + ct- ^ t 2 + •••)• 

(iv) /(A) = r (4±£) /r (£) with m(t) = ^ 3/2 - Here, 

~ y=t _1_1/2 (1 + ±ct - ±c 2 t 2 + •••)• 

(v) /(A) = T(aA + l)/T(aA + 1 - a). This lias density m(t) = T(1 _ a) ^*- t/a) i +a 
with asymptotics 

rtl (i\ - a 1+ “ /—I—a (-i , a-lf , 3a 2 -7a+2 + 2 , A 

myt) ~ pi-a) 1 ^ i + ~ t + —24P—t + j • 

All of the above densities are of exponential decay; we are not aware of densities 
that are of rapid but not exponential decay. 


Sample path properties. We next state the relationship between the asymp¬ 
totic expansion of the Levy density and sample path properties of the subordinate 
Brownian motion Bx t and the subordinator X t . 

A link with pathwise properties of the subordinate Brownian motion can be es¬ 
tablished via the Blumenthal-Getoor index that concerns the short-time behaviour 
of Bx t • This result is a simple consequence of mmm- 


Theorem 3.4. Let Bx t be the subordinate Brownian motion with subordinator X t 
whose Bernstein function satisfies Hypothesis \3.2l Then 


lim t~ 1/f3 sup\B Xs 

s<t 


with probability one. 


I 0 for all /3 > 2a 
oo for all (3 < 2a, 


The relationship with the subordinator is via Levy’s arcsine law. One can ex¬ 
press the order a and the coefficients po,pi, ... as expectations of suitable random 
variables: the asymptotics of m near t = 0 thus have a probabilistic representation. 

Theorem 3.5. Let the Bernstein function f satisfy Hypothesis 1 3. 2\ and let X t be 
the corresponding subordinator. For x > 0 define the first passage time strictly 
above x by T(x) = inf {t > 0|X t > x}. Then 

a = Ul ? + i E fe)' 

:r—>-0+ 

The lowest-order coefficient is given as 

Po = rpp 7 E (e“ A ' Yt ) 

for fixed t > 0 with similar expressions for the higher-order coefficients. 

Pseudodifferential operators and trace functional. We work in the algebra 
of SG- operators (sometimes known as scattering operators) as summarized in ]23j 
[24] , cf. also Section [2] for symbol spaces. 

We define the regularized trace functional TR exactly as in Section 2 of [24] with 
the exception that there is no integration with respect to x but only with respect 
to £, cf. equation m for an explicit expression. 
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A suitably shifted version of the generator A and the corresponding heat operator 
are classical pseudodifferential operators. 

Theorem 3.6. Suppose f is a Bernstein function satisfying Huvothesis \3. c A Define 
the operator A as in © and set A = A — m( 0, oo)/. Define coefficients 

(4) a k =-T(-ct +k)p k 

for k = 0,1,2,... Then the following holds: 

(i) The operator A is a classical elliptic pseudodifferential operator whose symbol 
<J(A) is in the class S c j*’ ( R") and has the asymptotic expansion 

OO 

(0~X>l£l 2a_2fc 

k =0 

in the sense of Definition ^. A 

(ii) The heat operator e~ tA is a pseudodifferential operator whose symbol cr(e tA ) 
belongs to sir 3 ' (R n ) with asymptotic expansion 

a(e tA )(£) ~ - [ai|e| 2 “ -2 + a 2 |£| 2 “~ 4 ] te~ ta °^\ 2a ± ■ ■ ■ 

in the sense of Definition \2.A 

The regularized zeta function. We define the regularized zeta function ((z) = 
TR where the complex powers of A are defined by functional calculus in the 

S'G-operators [23: . 

Theorem 3.7. Under the assumptions of Theorem } 3. 61 and with A = A— m(0, oo) I , 
the function (f(z) =TR is meromorphic on (D with at most simple poles at 

the points z k = {n — k)/2a for k = 0,1,2,.... The point z n = 0 is a removable 
singularity. 


Remark 3.8. The present exposition uses zeta functions as an intermediate step 
towards the heat trace. However, operator zeta functions are important in their 
own right in diverse branches of mathematics and physics, cf. m- We also mention 
for completeness that the residues of the regularized zeta function are expressed in 
terms of the noncommutative residue. This was originally defined by Wodzicki on 
closed manifolds extending of the work of Adler-Manin (cf. [33]) and plays a key 
role in noncommutative geometry. The residue depends only on the homogeneous 
component of order — n in the symbol expansion of A. 


(5) 


In the lowest orders, this residue becomes 
resjiz) = ^^al n/2a 
res C(“) = 0 

Z=Z\ 

res C(z) = — tT tt; sn u n 

z=z 2 K 2a (2-rr) n 0 

res C(z) = 0 


Qri —Z2~l 

a n a v z 2 


res C(*0 = h 


Wa -**~ 2 


Z = Z 4 

2 x"/ 2 


24 


24 ■ 


af — 2aoa2 

-Z2- 


2 a (2tt)" 2 u l u 0 

with fl n = 2 ( n '/ 2 ) surface area of the unit sphere in R” and the a k as in ©■ 
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The location and the residues of these poles are determined in terms of the as¬ 
ymptotic expansion of m near t = 0. By Theorems 13.41 and 18.51 we can express this 
information probabilistically in terms of the subordinator. 


Heat trace expansion. The pole structure of the regularized zeta-function 
determines the short-time asymptotics of the generalized heat trace TR(e~ tA ). 

Theorem 3.9. Under the assumptions of Theorem, 1,9. 61 and with A = A—m( 0, oo) I, 
the asymptotics as t —> 0 + of the generalized heat trace TR(e ~ tA ) are given as 
follows. 

(i) If a is rational, there are constants c k and c k such that 

OO OO 

TR(e~ tA ) ~ Ckt~ {n ~ k)/2a - ^ log t. 

fc=0 k=l 

(ii) If a is irrational, we have 

OO 

TR{e~ tA ) ~J2ckt- {n - k)/2a 

k—0 

where c k = T(zz±) res C(*)- 

z=(n—k)/2a 


Since A is merely A shifted by a constant we note the immediate 


Corollary 3.10. Under the assumptions of Theorem, 1. 9. .91 the heat trace expansion 
of A is given as follows. 

(i) a rational: there are constants c k and c k such that 


OO OO 

TR(e~ tA ) ~ e ~™(o,oo)t ^ Ck t- { - n - k) /‘ 2a - e"™ (0lOo)t ^ c k t k log t. 

k =o fc=1 

(ii) a irrational: we have 

OO 

TR(e~ tA ) ~ e -™(°.o°)t^ Cfet -(n-fc)/2a 

fc=0 


for constants c k . 


Note the strikingly different behaviour for a rational and irrational with the 
appearance of logarithmic terms. For dimension n > 2 we explicitly give the lowest- 
order terms (the case n = 2 is illustrated in Section 0: 


_ T(n/2a) —n/2a 


2a (27r j n a 0 


( 6 ) 


Co = 

ci = 0 


_ r((n—2)/2a) n n 

02 2 a 


-(n-2)/2a-l n-2 


' (2ti-)" u 0 ai ' 2 a ’) 

where the a k are as in 0). The lowest-order term co also appears in [7] and [5B], cf. 
Remark IQ The lowest-order non-logarithmic terms of the asymptotic expansion 
read 


(271-)" 


1 

2a 


n%)a- n/ 2 a t- n / 2a - r (^) a - ( «- 2)/2 “-i ai ^i 


—2j-—(n—2)/2o 


geometry 
of lR n 


dimension n and probabilistic properties of Bx t 
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Note that the terms in the heat trace expansion combine the two main aspects 
of our subordinated Brownian motion: the flat Euclidean geometry in which the 
process lives and information about the Levy density. 

Remark 3.11. We briefly comment on extensions and limitations of our approach. 

(i) The approach works whenever the generator A (or a suitably shifted genera¬ 
tor) belongs to an algebra of classical pseudodifferential operators that has a 
trace functional and contains the complex powers A~ z and the heat operator 
e~ At . Thus, on a closed manifold M the same calculations can be carried 
out demonstrating how both the classical heat invariants and certain prob¬ 
abilistic information of the subordinator appear explicitly in the heat trace 
asymptotics. 

(ii) The definition of the regularized trace TR. and the meromorphic extendibility 
of C(z) require classical symbols. This limits the class of Markov processes 
that can be treated by this method. 

(iii) Since the symbol calculus of pseudodifferential operators works only modulo 
smoothing operators (symbols of order — oo), the asymptotic expansion of 
a(A ) cannot see the behaviour of the Levy density m{t) for values of t away 
from 0. Thus, one would have to work in a different operator algebra if one 
wanted to capture the behaviour of m at both t = 0 and t = oo. 

The present exposition naturally leads to further questions that are beyond the 
scope of the present exposition and are the subject of further research. 

(i) Is there a probabilistic characterization, e.g. in terms of sample path proper¬ 
ties, of our class of Bernstein functions? 

(ii) How can one recover the ’’long end” of the asymptotics of to, i.e. as t —> oo 
in the spectrum of A? 

(iii) What is the probabilistic significance of the logarithmic terms in the heat 
trace and the dichotomy rational/irrational a? 

4. Proofs of the key results 

This section contains the building blocks needed for the proof of the key results. 
We have structured the arguments into the following subsections. 

(1) Growth, regularity and asymptotics of a class of Bernstein functions 

(2) Sample path properties 

(3) Constructing a classical pseudodifferential operator 

(4) The zeta function and heat trace of the shifted operator 

4.1. Growth, regularity and asymptotics of a class of Bernstein functions. 

We derive lower and upper bounds for our class of Bernstein functions, establish 
their smoothness and give explicit asymptotics. 

As a general remark, note that by Theorem 3.9.29 of [18], any Bernstein function 
can be continuously extended to [0, oo) and we tacitly assume that all our functions 
are defined on this larger interval. 

The following proposition suggests that we can construct a pseudodifferential 
operator with elliptic symbol given in terms of / with an asymptotic expansion 
into homogeneous terms. This will be taken up in detail in the next subsection. 

Proposition 4.1. Let f be a Bernstein function satisfying Hypothesis \S.2\ Then 
the following holds. 
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(i) We have f £ C' oo ([0, oo)), i.e. f is smooth at the origin, 
(ii) The function f has an asymptotic expansion 


(7) 


/(A) ~ m(0, oo) - ^2 T(-a + k)pu A 1 


a—k 


k—0 


as A -A oo where m( 0, oo) = / 0 °° (m(t) — pot 1 “) dt. 
(in) The derivative has an asymptotic expansion 


( 8 ) 


/ { 0 (A)~ (-l) i+ 1 ^r(-a + « + fc)p fc A c 


— l — k 


k—0 


as X ^ oo for any l = 1,2,... 

There is an R > 0 swc/i that 

-±T(-a) Po X a < /(A) 

/or X> R. 

(v) There are constants Ci with |/^(A)| < Ci\ a ~ l for A > 1 and l = 0,1, 2,... 

Remark 4.2. If we allow more general Bernstein functions of the form /(A) = 
a + bX + J 0 (1 — e~ xt )m(t)dt with a > 0 and b > 0, then the expansion in (ii) has 
the additional terms a + bX. Also, the lower bound in (iii) reads i bX. 

The key argument for the proof of assertions JnJ and (lull) , from which (IrvT) and 
(jvj) follow immediately, is a version of Watson’s Lemma that we state without proof. 

Proposition 4.3 ([10], Section 4.1). Let g : (0,oo) — > 1R be a bounded and locally 
integrable function such that as t —> 0 + one has 


9(t ) ~ 5>ki a ». 


k—0 


Assume pk G ak > — 1 for all k and ak —> oo monotonically as k oo. Then 


poo 00 

/ e~ M g(t)dt ~ r(l + a k )p k X 

"'° k—0 


— ak — 1 


as X 


00 . 


Proof of Proposition E3 0 The smoothness of / is a consequence of the rapid 
decay of m as t —> oo since then / 0 °° t l m(t)dt < oo for l = 1,2,.... By Theorem 
3.9.23 of [TS] the limit /(0) = lim t _ >0 + f(t) exists so that / is continuous at the 
origin. A direct calculation interchanging differentiation and integration shows that 
/ is differentiable on the half-line [0,oo). Moreover, any derivative of / on (0, oo) 
(given by differentiating the Laplace transform of m) can be continuously extended 
to 0 as all higher moments of m exist due to its rapid decay. We can repeat these 
arguments for any derivative of / so that / is smooth on [0,oo). 

JnJ This follows from Watson’s Lemma. We consider 


poo poo poo 

/ (1 -e~ xt )m(t)dt= (l - e~ xt )p 0 t~ 1 ~ a dt+ (l - e~ xt ) ’ 

Jo Jo Jo 


l(t)dt 


with m(t) = m(t) — pot 1 a . The first integral can be evaluated and yields 


( 9 ) 


poo 

/ (l - e~ xt ) p 0 t~ 1 ~ a dt = r(1 ~ a) PoA a = -T(-a)p 0 A“. 
JO 
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The second integral can be expanded using Watson’s lemma where we have g(t) = 
fn(i), ao = 0 and a k = — 1 — a + k for k > 1. This yields 


( 10 ) 


p OO pOO pOO 

I (l — e~ xt ) m(t)dt = / m(t)dt — I e~ xt m(t)dt 

Jo Jo Jo 

oo 

~ m( 0, oo) — ^ T(-a + k)pk A“ 


—k 


fc= 1 


for A —> oo. Adding © and m proves the assertion. 
m Note that any derivative of / is given in the form 

POO 

/ W (A) = (-1)' +1 / e~ xt t l m{t)dt. 

J o 

By Hypothesis 13.21 the map t l m(t) is integrable for any l > 1: at t = 0 it grows 
like and at infinity it is of rapid decay. So by Watson’s Lemma 


/ W (A) ~ (-1)' +1 ^r (-a + l + k)p k X a ~ l ~ k 

fe=o 


for A —> oo. 

(IrvT) The bound follows from the lowest-order asymptotics just derived. We find 


lim A “/(A) = —T(—a)p 0 - 

A—>-oo 


This means that for every e > 0 there is some R > 0 such that |A “/(A) — po\ < e 
for all A > R. Now set e = — T(— a)po/2 so that 

(11) -|r(-a)po<A-“/(A)<-|r(-a)po 

for all A > R. This proves the claim. 

(jvj) The bound for l = 0 follows from (fill) . The claim for l > 1 follows by 
considering lim^oo X~ a+l (X) = (—1)* +1 T( —a + l)po by (Iml) and arguing as in 
(1A1) . One can have R = 1 by choosing Ci sufficiently large. 

□ 


4.2. Sample path properties. We prove the assertion on the growth of the sam¬ 
ple paths of Bx t ■ 


Proof of Theorem \S.4\ This is a simple application of the results in 2S]. Denote 
by = /(|£| 2 ) the symbol of the process Bx t - As in Example 5.5 of [2H] we 
compute the Blumenthal-Getoor-type indices 


/Soo = inf |r > 0 
Soo = inf •( r > 0 


lim = 0 

Ifl-yoo Kl 


lim inf = 0 l. 


l£|->oo 


l€l 


By Proposition 14.11 |u| , / has the asymptotic expansion 


OO 

(12) /(A) ~ m( 0, oo) -Y, r (-« + k) Pk X a ~ k 

k -0 

as A -I oo. Thus Poo = Soo = 2a and the claim follows from Theorem 4.6 of 
[28]. □ 
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We also give a proof of Theorem 13.51 yielding a probabilistic interpretation of 
the coefficient cio in the asymptotic expansion of m under Hypothesis 13.21 This 
interpretation will hinge on Levy’s arcsine law. To this end recall the notion of 
regularly varying functions. 


Definition 4.4 ([9], Chapter 0.7). A measurable function / : (0,oo) —> (0, oo) is 
regularly varying at oo if for every r > 0 the ratio /(tA)//(A) converges in (0, oo) 
as A —> oo. Indeed, there is a real number p, called the index such that 


lim 

A—>-oo 


fW 


= T P 


for every r > 0. 


Proof of Theorem 1 8. -51 By equation m we find for r > 0 that 


(13) 


f(T A) 

fW 


—> T 


a 


as A -A oo since in m the highest power in A is of the order a > 0. This means 
that / is regularly varying at oo with index a. 

By Levy’s arcsine law (Theorem III.3.6 of [5]) we have 


(14) 




where T(x) = inf {t > 0\X t > x} is the first passage time strictly above x. 

The assertion that the coefficient po can be described probabilistically is clear 
since A”“/(A) —> —T(—a)po for A —> oo. As both ao and / can be described 
probabilistically in terms of X t via (fl4l) and E ( e ~ XXt ) = e _t ^ (A ), respectively, so 
can po: 

~T(-a)po= lim A~ Q /(A) = lim -A" Q ± logE (e“ A ' Yt ) 

A—^oo A —y oo 

for t > 0. Likewise for the other coefficients pk■ □ 


4.3. Constructing a classical pseudodifferential operator. We shift the the 
generator A by a suitable constant to obtain a classical pseudodifferential operator 
A = A — rh(0, oo)/. The classicality is crucial in order to define the regularized zeta 
function and generalized heat trace. The heat kernel will be a classical pseudodif¬ 
ferential operator if A is invertible and parameter elliptic with respect to a sector 
in the complex plane, cf. Definition 12.31 (fiH) . 


4.3.1. Classicality of the symbol. It is easy to see that a suitably adjusted version 
of /(l£| 2 ) is a classical symbol. 

Proposition 4.5. Assume Hvvothesis \3.2\ and define a function 

ct(£) = /(ICI 2 ) — TtT(0, oo). 

Then d is a symbol in 5 2a,0 (R n ). Moreover, it belongs to S^’°( 1R"), i.e. it is 
classical with asymptotic expansion 

OO 

(15) m~-Y,n-a+k)pm 2{ - a - k) 

k—0 


in the sense of Definition ] 2. SI 
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Proof. 1. We first show that a £ 5 2 “ , 0 (IR"). By Definition 12.11 it suffices to show 
that /(|£| 2 ) is a smooth function on R" and that for all multi-indices /3 £ Ng there 
is a constant Cp with 


(16) 


^/(|£| 2 ) 


< Cp(0 


2a-|/3| 


for any £ £ 1". Recall that (£) = y/1 + |£| 2 . 

Now, /(|£| 2 ) is smooth as / is smooth on [0,oo) by Proposition 14.II (HI). 

To obtain the estimates in (fl6l) recall Faa di Bruno’s formula for the chain rule 
in the version of equation ( 6 . 8 ) of [32]. Let /3 £ N" be a multi-index. Then 

3f/(|£| 2 )= E c p , 7l ,... i7p (d^f) (|£| 2 )c^ 1 |£| 2 • • • <9j p |£| 2 , 

7i H-I-7 p=/3 


where the sum runs over all partitions of the multi-index /3 into sums of nonzero 
multi-indices 71 ,..., for p = 1,2,... As |£| 2 is a symbol in S , 2 , 0 (R n ) we have 


«| 2 <c 74 (£) 2 - |7i| 


for some C 7i . Also, | 7 i| H-+ | 7 P | = |/3| by construction so that overall 


«| 2 -'-tt| 2 


— C7i,...,7p 


(0 


2p-|/3| 


for a constant C 71j .., )7 . From Proposition 14. II (frvT) we find 


\(dif)(\e)\<c P \e {a ~ p) <c P (o 2[a ~ p) 

for some constants C p so that overall m follows. 

2. Now prove the asymptotic expansion (fl5l) . Let x G C'°°(R”) be a cutoff 
function which is equal to 1 outside |£| > 1 and equal to zero for |£| < 1/2. We 
must show 


N-l 

<7(0 -m(0,oo)+ £ xm -a + k)p k \^ a -V £ S 2a ~ N ’°OR, n ) 

k —0 


for any N > 1. 

Suppose |£| > 1 and consider /(|£| 2 ) — Tfi(0, 00 ) + J2k=o r (— a + k)Pfc|£| 2 ^ a- ^- 
We can rewrite this as /at(|£| 2 ) with 


POO / N-l \ 

In( A) = (1 - e~ xt ) ( m(i) - 53 J dt. 

Now repeat the arguments of Proposition 14.11 (IrvT) to see that 

fjp( a ) <c hN x a ~ N ~ l 


for some constants Ci , n- An argument as in 1. shows that Jn belongs to S 2a ^(R") 
proving the claim. □ 
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4.3.2. Ellipticity of the symbol. We note the (parameter-)ellipticity of a shifted 
version of the Bernstein function. 

Proposition 4.6. Assume Hypothesis 1 8. 2\ and set cr(£) = /(|£| 2 ) — m(0,oo). Then 
a is elliptic. Moreover, for 9 £ (7r/4, n/2) the symbol a is A-elliptic. 

Proof. 1. For the ellipticity we must show that there is a C > 0 with 

(17) C(0 2a < <7(0 

for any £ e R n . 

By Proposition 14.11 (IrvT) there is an R > 0 such that for |£| 2 > R we have 
Ci|£| 2a < <7(0 where C\ = — |T(— a)po- Without loss of generality R > 1 so that 
for |£| 2 > R we have C 2 (£) 2a < <r(£) with C 2 = (i)“Ci as ±(£) 2 < |£| 2 for |£| > 1. 
There is an R' > R such that 

IC 2 (0 2a <C 2 (0 2 “-m(0,oo) 

< /(ICI 2 ) - m(0,oo) 

= o -(0 

for |^| 2 > R'. 

If |^| 2 < R' we have 

—m(0, oo)(i?') _2 “(0 2a < — m(0, 00 ) < cr(0- 
Thus, (fl7l) holds with C = m.\n.{^C 2 \ — m(0, oo)(i?') _2a }. 

2. Now for the parameter-ellipticity. We must show that 

l(A-<r(0r 1 |<C<0" 2a 

for some C > 0 and all £ £ IR". Equivalently C' _1 (£) 4 “ < |A — cr(£)| 2 . Let A £ A 
with A = Ai + i\ 2 and Ai £ IR. We distinguish cases according to the sign of Ai. 
Case Ai < 0: we find 

|A - a{x, £)| 2 = (Ai - <j (£)) 2 + A 2 

> ^(0 2 

> c'(0 4a 

for some constant C' by ellipticity and since Ai and —a have the same sign. 

Case Ai > 0: here A 2 > Ai since we assumed 9 > 7r/4. Thus, 

|A — d (£)| 2 = (Ai — cr (£)) 2 + A§ 

> (Ai - d (£)) 2 + A 2 

> ^(0 2 

> \C'{£) 4a 

for some C' by ellipticity. Here we use the fact that the function g{x) = (x — c) 2 +x 2 
for c > 0 has a minimum located at x = c/2 with value g(c/ 2) = c 2 /2. □ 

4.3.3. The spectrum of the operator. We collect useful facts about the spectrum of 
the operator A = A — m(0, 00 ). 

Proposition 4.7. Assume Hvvothesis \3.2\ and set <r(£) = /(|£| 2 ) — m(0,oo). Let 
A be the pseudodifferential operator with symbol a. 

(i) A is essentially selfadjoint, i.e. has self-adjoint closure A. 
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(ii) The spectrum of A is contained in the positive real line so that in particular 
A is invertible: Sp C [—m(0, oo), oo). 

Note that the closure of A generates an analytic semigroup which by abuse of 
notation we denote by e~ tA . The heat operator is a classical pseudodifferential 
operator. 

Corollary 4.8. Under Hypothesis 1 3. 2\ the operator —A extends to the generator of 
an analytic semigroup e~ tA on L 2 (R”). The heat operator e~ tA is a classical pseu¬ 
dodifferential operator with symbol u(e~ tA )(f) = e~ ta ^^ 1 belonging to <S'd oo,0 (lR”). 

Proof. The spectrum of A is entirely contained in the interval [—m(0, oo), oo) by 
Proposition 14.71 JnJ). Thus, the spectrum of —A is contained in a sector in the left 
half of the complex plane and the analytic semigroup e~ tA exists by Theorem II.4.6 
of j Tfij . 

The fact that e~ tA is a pseudodifferential operator follows from a simple appli¬ 
cation of Theorem 4.1 of [24] whose hypotheses are satisfied: the symbol d is of 
positive order in £ and order 0 in x. Moreover, it is A-elliptic by Proposition 14.61 
and the resolvent of A exists in the whole sector A as 0 ^ Sp(A) by Proposition 14.71 
©. ' □ 

Proof of Proposition \ f. 7 \ (i) By standard arguments (cf. for example Chapter 1 of 
[30] or Chapter 2.1 of .12]), A is closable. Also, A is symmetric, i.e. (Au,v) = 
( u,Av ) in L 2 (R n ) for any u, v £ <S(R”) since d is real-valued and independent of 
x. By Proposition 1.3 of [12], A is selfadjoint if and only if A is symmetric. 

(ii) Spectrum of A: clearly, Sp ^A^ C [0,oo). Proposition 3.10 (ii) of [30] implies 

that for a selfadjoint densely defined operator T we have A ^ Sp(T) if and only if 
||(T — A)u|| > c_\||u|| for some constant c\ > 0 and all u in the domain of T. In our 
case note that for u £ <S(R") we have ||Au|| 2 > (m(0,oo)) 2 ||u|| 2 by Plancherel’s 
theorem as the symbol of A is bounded below by —m(0, oo) so the claim follows. 

□ 

4.4. Zeta function and heat trace of the shifted generator. For a classi¬ 
cal pseudodifferential operator we compute the regularized zeta-function and the 
generalized heat trace. 

Proof of Theorem 1 3.6 1 (i) The assertions follow from Propositions 14.51 and 14.61 
(ii) We set ak = —pkU(—a + k) for k = 0,1,2,... for short and apply Theorem 
4.1 of [24]. Let A be a sector as in Definition 12.31 (fn|) with 0 as in Proposition 14.61 
To find the symbol of the heat operator we first construct the symbol expansion of 
the parameter-dependent parametrix using the algorithm of Proposition 3.2 of [231 
for symbols independent of x. Recall that by a parameter-dependent parametrix 
of a pseudodifferential operator A with symbol <r(£) £ S r '° we mean an inverse 
modulo smoothing operators, i.e. operators with symbol in iS' _oo ’ 0 (R n ). If we 
denote the symbol of the parametrix by 6(£, A) then on the level of symbols we 
require that (A — cr(£))b(£, A) — 1 £ S -JV,0 (IR, n ) for any IV £ IN. The ansatz b(£, A) ~ 
b- r (f, A) + &_r-_i(£, X) + b- r - 2 (f, A) -)-leads to the relation of formal power series 

[(A — <5y(£)) + <5y_i(£) + • • • ] [&-r(£, A) + b- r - 1 (£, A) + • • • ] = 1. 
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Roughly speaking one collects terms according to the degree of homogeneity in £ 
where A — ay is treated as order r. 

In our case with r = 2a the first terms are given as 

b-2a(£, A) = 
b~ 2c*-i(£, A) = 0 

(18) b_2oi—2 (£: A) = ( A _ 1 ( lfJ| ? |2Q ! )2 > 

6-2a-3(£i A) = 0 

L ^ M ^kl 2 “- 4 . «?l€l 4 “- 4 

0_2a-4(.t, A; - ( A -aokP») 2 + (A-a 0 k| 2 “P 

The terms of order — 2a— k with k odd all vanish because the asymptotic expansion 
of <7 contains no terms of order 2a — k with k odd. 

We can then determine the symbol expansion for the heat operator using the 
above parametrix to obtain 


—tA 


)«) ~ E 


fc=o 


2m 


' dA 


e~ tX b-2a-k{^,, X)d\ 


by Theorem 4.1 of [!M], where dA is a parametrization of the boundary of the sector 
A. To the above lowest order terms, the symbol is thus given as 

<j(e tA )(Z) ~ e - taol « |2 “ + [—ai |£| 2 “" 2 -a 2 |e| 2 “" 4 ] te~ ta °^ a 

(19) +\u\\Z\ 4a - 4 t 2 e" too|?|2 “+••• 

which follows by integration by parts (detailed computations in Section [5]). □ 


We can perform a similar computation for the ((-function of the operator. 


Proof of Theorem \3.7\ Let A be a sector as in Definition ^. 31 (ln| with 9 as in Propo¬ 
sition m For the construction of the complex powers of A with classical symbol 
a we must have that 

(i) A — A is invertible for all 0 ^ A £ A and 

(ii) A = 0 is at most an isolated spectral point. 

These assumptions are satisfied in our case by Proposition l4.7l([iI)) . Then define the 
complex powers A~ z for 2 £ <D, Re z > 0 by a Dunford integral 


A~ z 



A -2 (A — A)~ 1 d\ 


for Re z > 0 where A e = AU{|A| < e} and dA e is a parametrization of the boundary 
of A e with the circular part traversed clockwise. We refer to [23], Section 3.2 to see 
that A~ z is a pseudodifferential operator in S' 2 “ 2: ’ 0 (R n ) with symbol <t(£; 2 ) and 


a (& z )=2h [ x Z ( X - (J (0) 1 d\■ 

Here, = {2 £ (D \ A|i(£) r < | 2 | < c(£) r } for a suitable constant c > 1 such 
that ct(£) £ for any £ 6 R", cf. Lemma 3.2 of [23]. The symbol is defined 
modulo smoothing symbols in S' -oo,0 (R ra ) depending analytically on 2 . 

The asymptotic symbol expansion of A~ z is given for Re 2 > 0 as 

a (^ _Z ) ® ~ z) + CT- 2 a^- 2 (C; z) + Cr_ 2a z-4(£; z) H- 
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with 

<T-2a.(^) = ao Z |r 2 “* ] 

(20) a. 2az - 2 (t;;z) = -ao Z - 1 a 1 z\{\- 2az - 2 l 

a- 2a z- 4 {.£,]z) = [-aQ Z ~ 1 a 2 z+ \o.q Z ~ 2 a\z(z + 1 )] \£\- 2az ~ i ) 

obtained by integration by parts (exemplary computations in Section O . 

Now for the regularized zeta function ((z) = TR(A~ Z ). By Theorem 2.7 of [24] . 
this function is analytic for Re 2 > nj2a and can be meromorphically extended to 
the entire complex plane with at most simple poles at the points Zk = (n — k) /2 a 
for k = 0,1,2,.... The behaviour of £(z) near Zk is described as 

( 21 ) 2 a(z-z k ) n [ cr -2az+2otz k -n (?! z)du (£) 

with w the surface measure on the unit sphere in R". Compared to (2.10) in [24] all 
terms with averaging in the ^-direction are omitted, the finite-part integral vanishes 
anyway on operators with symbols independent of x. 

□ 

Based on the singularity structure of the regularized zeta-function we can give 
the short-term asymptotics of the generalized heat trace. 

Proof of Theorem \3.9\ This is an application of Theorem 4.3 in [24]. It relies on a 
mathematical folklore written up in Section 5 of Hzj. The singularity structure of 
the generalized heat trace is related via the Mellin transform to r(z)£(z) so that 
we need to consider the pole structure of T(,2)£(,3:). Indeed, for 2 a rational 

oo oo 

( 22 ) T(z)C(z) ~ z-(n-k)/2a + [z-(-i)] a 

k—0 1=1 

with coefficients Ck and q as in the assertion. We have double poles where the poles 
of T(z) and £(*) coincide which can only occur for (n — k)/2a a negative integer 
denoted by l in equation (]22l) . 

For 2a irrational, r(z)C(z) has only simple poles T(. 2 )C(z) ~ z _( n c \)/ 2a 

with Ck as given. □ 

5. Worked example: the relativistic stable Levy process 

To illustrate the approach and to show that it is analytically tractable we consider 
the Bernstein function /(A) = y/X + 1 — 1 in dimensions n = 2 and n = 3. The 
generator of this process is usually denoted by — H where H = (—A + 1) 1//2 — 1 is 
the relativistic Hamiltonian with zero potential and mass 1 and A is the Laplace 
operator. The behaviour of the regularized zeta function and the generalized heat 
trace is different in these dimensions with the appearance of logarithmic terms in 
dimension 3. 

5.1. Dimension-independent considerations. We first derive several proper¬ 
ties of the Bernstein function and the associated pseudodifferential operator that 
are independent of the dimension n. 
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Asymptotic expansion of the Bernstein function. The Levy density of 
the Bernstein function is given by m(t) = as L(l/2) = ypK. From the 

Taylor series for e~* near t = 0we find 

m (t) ~ 27¥ i_3/2 (! - t + it 2 - ± ■ • • ) . 

1 ( _ 1^1 

In our previous notation this yields a = —^ and pk = v k\ f° r ^ = 0,1, 2,... 

We apply Proposition l4.ll © to find the asymptotics of /. With m(t) = m(t) — 
^A=f -3 / 2 we find m(0, oo) = — 1. Thus, 

/(A) ~ A 1 / 2 - 1 + iA" 1/2 - |A- 3/2 + ^A" 5/2 T • • • 
which agrees with the Taylor expansion of / in 1/A. 


The associated pseudodifferential operator. The associated pseudodiffer¬ 
ential operator A has the symbol cr(£) = (|£| 2 + l) 1 / 2 — 1. We consider the shifted 
operator A = A + I with symbol 

*(o = n\e) +1 ~ i£i + nr 1 - m ~ 3 + ^ir 5 t •• • 

which is a classical symbol in 5' 3 ; ’°( 1R 2 ). 

Let A be the sector A = {z £ (D|0 < arg(z) <2n — 9} for a 9 € (7r/4,7 t/ 2). A 
quick calculation shows that the symbol is A-elliptic, i.e. that <r(£) does not take 
values in A and ^=(£) < |A — cr(£)| for all A G A, £ G R 2 . 


Parameter-dependent parametrix. For the parameter-dependent parametrix 
we apply the algorithm of the proof of Theorem 13.61 and find 

1 


b- i(&A) = 


6_a(^,A) = 2 
A) = 


A-iei 

Hr 1 


(A-|?|) 2 

(Hr 1 ) 2 


iltl-3 


m 


(A^kl) 3 (A — |£|) 2 ’ 
which agrees with m- Putting everything together we arrive at 


b(Z, A) 


1 




A- ICI (A-|el) 2 

for the parameter-dependent parametrix. 




1 ICI-3 


m 


(A — |£|) 3 (A-I^l) 2 


+ 


Symbol expansion of the heat operator. The first few terms in the asymp¬ 
totic expansion of the symbol of the heat operator are given by 

a (e- tA ) (0 ~ e-^l + + A|£|- 3 ] te~™ + fir 2 ^"* 1 * 1 + • • • 

in the lowest orders. 


The regularized zeta function. The lowest-order terms in the asymptotic 
expansion of the complex powers is given by 


(23) <t (a-*) (o ~ ier - *- 2 +§(* 2 +2*)ir * -4 + • • • 
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for the complex powers. 

5.2. Dimension n = 2. The regularized zeta function has at most simple poles at 
the points Zk = 2 — k for k = 0 , 1 , 2, ... and the residues are given in terms of the 
homogeneous symbol of order — n = —2. Let u> be the surface measure on the unit 
sphere in R 2 . 

k = 0 : near zo = 2 — 0 = 2 we find that 

cw = / ur’-MO = A - A 

k = 1 : near z\ = 2 — 1 = 1 we see that £(z) is analytic as the asymptotic expansion 
(1231) contains no homogeneous term of order —2. 
k — 2 : near z 2 = 2—2 = 0 , the homogeneous term of order 2 is given by — ^z \^\~ z ~ 2 
which vanishes at Z 2 = 0 so that £(z) is analytic there. 
k = 3: near z 3 = 2 — 3 = —1 we find that £( 2 ) is analytic as the asymptotic 
expansion (1231) contains no homogeneous term of order —2. 
k = 4: near Z 4 = 2 — 4 = —2 we have 

t( z ) = ^+rk [ I [^ 2 + 2 z] 2= _ 2 l£r 2 <M£) = °, 

so that £(z) is analytic there. 

Hence to lowest orders the only pole of £(z) is at z = 2 with residue -A. 


Heat trace expansion. We finally compute the coefficients of the heat trace 
expansion. The ansatz is to write 




00 


E 


Cfc _I_ <-k 

z—(2—k) ^ (z-(2 -k)) 2 


+ g{z) 


where g is an entire function which is of no concern to us. The T-function has only 
simple poles located at Zk = —k for k = 0,1,2 ,... with residue res 2 = _/ c T(z) = 
- . An explicit computation yields the following: 

c° = resT(2)C(^) = T(2) • ^ = ^ 

ci = 0 

C2 = 1 • C(0) 

c 3 = res T{z)C(z) = . ((- 1 ) = -((- 1 ) 

z=—l 

c 4 = res T (z)C(z) = klL . ((- 2 ) = |C(~2) 

z=—2 


with Co = Ci = £2 = £3 = C 4 = 0 as there are no double poles of r(z)£( 2 i) at these 
points. 

Hence the desired heat trace expansion for A reads 

TR (e~ tA ) ~ ±t~ 2 + C(0) - C(-l )^ 1 + K(-2 )* 2 + • • • 

in agreement with ([ 6 ]). Note that the heat trace expansion of the relativistic Levy 
process TR ( e ~ tA ) can be obtained via the relation TR ( e ~ tA ) = e t TR as 
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A = A + / so that 

TR (e~ tA ) ~ + C(0)e 4 - C(-l)eV + iC(-2)e‘t 2 + ■ ■ ■ . 

In dimension n = 2, the zeta-function is regular at —k for k = 0,1,2,... (and so 
there are no logarithmic terms in the heat trace expansion) as n = 1 /a, i.e. there 
is a algebraic relation between the dimension n and the parameter a. 

5.3. Dimension n = 3. The regularized zeta function has at most simple poles at 
the points Zk = 3 — k for k = 0,1,2,... with the following residues. Let ui be the 
surface measure on the unit sphere in 1R 3 . 

k = 0: near zo = 3 — 0 = 3 we find that 

C(*) = • 3 / I Z\~ 3 M0 = 6^ ■ 1=3 

k = 1: near z\ = 3 — 1 = 2, the function Q{z) is analytic. 

k = 2: near Z 2 = 3 — 2 = 1, we obtain 

C(*) = • I / (-3) • 1 • ir 3 <M£) = 

•' 141=1 

fc = 3: near Z 3 = 3 — 3 = 0, the function ((z) is analytic. 

k = 4: near Z 4 = 3 — 4 = — lwe find that 

CM = ■ I [ I [(-1) 2 + M-1)] = 

Heat trace expansion. Explicit computations as above yield the following 

c ° = reg r (z)C(ar) = T(3) ■ ^ ^ 

ci = 0 

c 2 = rraT(«)C(«) = r(l) • = - T 

c 3 = res T( 2 ;)C(z) = 1 • C(0) = C(0) 

z=0 

C 4 = 0 . 

Also, co = ci = c 2 = 63 = 0 and due to the double pole of r(z)()(z) at Z 4 = —1 

C4 = ^ • (-3^r) = 

The desired heat trace expansion for A reads 

TR (e~ tA ) ~ 3 ^eH ~ 3 - j^eH - 1 + C(0)e‘ - j^eH log* + • • • , 
for the relativistic stable Levy process in R 3 . 

Remark 5.1. We remark that [7] and [26] consider a related case when an a-stable 
relativistic process moves in a bounded region D in Euclidean space of dimension 
n. The authors assume different degrees of regularity of the boundary dD. In our 
notation, the authors consider the subordinator /(A) = (A + cV a ) a — c. 

For c = 1 the two lowest orders of the heat trace are given in 13 by 

(24) vol(D)Ci {t)e t t~ n ^ 2a - C 2 (t)vol(0D), 
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and in [26] by 


(25) 


vo 1(D) Ci 



C' 2 vo\(dD)t- in - 1 '> /2a , 


where l is the largest integer less than 1/a. Ein either case, vol(D) the volume of 
D with respect to Lebesgue measure and vol(<9-D) the surface area. 

In (l24l) . C\ and C% are functions such that 


(26) 


C\ (t) —» 


n n F(n/2a) 

(2ir) n 2a 


as t -» 0 and C^it) is bounded in the order of e 2t t 1 l/ 2a . Both C\ and Ci 
can be expressed in terms of probability densities and first exit times, respectively. 
In (l25l) . C'i is a constant given by C[ = ■ The constant C' 2 reflects the 

boundary geometry of D. 

Note that the limit in (l26l) (and hence C[) is exactly the coefficient Co in our 
heat trace expansion ([6]) for the given subordinator. 
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